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This paper presents a collection of formulas that are necessary for the 
treatment of radiation and mode conversion phenomena of imperfect 
asymmetric slab waveguides. The coupled mode theory of dielectric wave- 
guides is briefly reviewed, and general expressions for the coupling coeffi- 
cients are given. The field expression of the guided and the radiation TE 
and TM modes of the asymmetric slab waveguide are stated, and are used 
to derive formulas for the coupling coefficient for slight core boundary 
irregularities. 

I. INTRODUCTION 

Mode coupling phenomena and radiation losses caused by core- 
cladding interface irregularities have been studied extensively for 
symmetric slab waveguides and for round optical fibers. 1-7 These 
results are not immediately applicable to the asymmetric slab wave- 
guides used in integrated optics circuits. It is the purpose of this paper 
to collect the formulas for the normal modes of the asymmetric slab 
waveguides, and for the coupling coefficients between guided modes 
and guided and radiation modes caused by core boundary irregularities 
of these waveguides. 

The coupling coefficients between guided modes are useful for the 
design of distributed feedback sections for lasers and for an evaluation 
of unintentional mode coupling caused by core boundary roughness. 
The results collected in this paper are further necessary for the evalua- 
tion of radiation losses caused by core boundary roughness. 

Because of the many parameters that enter into the theory, it is 
impossible to evaluate the formulas in graphical form for all cases of 
practical interest. This paper is thus a collection of the required 
formulas which the reader can use to evaluate his particular problems. 
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II. SUMMARY OF THE COUPLED MODE THEORY 

The coupling theory is based on an expansion of the solution of 
Maxwell's equations in terms of normal modes. The general theory of 
the mode expansion and the derivation of the coupling coefficients have 
been published by A. W. Snyder. 2,3 His theory is based on local normal 
modes. Local normal modes resemble the modes of the ideal asymmetric 
slab waveguide with perfect core boundary. However, the boundary 
of the perfect guide is allowed to change in such a way that it coincides 
with the actual deformed core boundary at the particular point z along 
the waveguide axis at which the coupling coefficients are to be eval- 
uated. Since the waveguide width parameter is no longer a constant, 
the local normal modes are not themselves solutions of Maxwell's 
equations. They must be superimposed with z-dependent expansion 
coefficients to form such a solution. The fact that these modes form a 
complete orthogonal set and coincide with the modes of a fictitious 
waveguide, the width of which is locally (at the point z under consider- 
ation) the same as that of the deformed waveguide, explains the name 
"local normal modes." It is also possible to express the general field 
in terms of the modes of the ideal waveguide, the constant width of 
which differs slightly from that of the actual waveguide. This expansion 
suffers from convergence difficulties that are caused by the fact that 
the normal components of the electric field are discontinuous at the 
core boundary. The modes of the ideal guide are discontinuous at the 
dielectric interface of the ideal guide which does not coincide with that 
of the actual guide. The expansion in terms of ideal modes of the 
waveguide is thus discontinuous term by term at a point where the 
entire series must be continuous, and furthermore, it must describe a 
discontinuous field at the interface of the actual waveguide at a point 
where each individual term of the expansion remains continuous. The 
expansion in terms of local normal modes, on the other hand, describes 
the field discontinuity with a series, the individual terms of which 
are discontinuous in just the right way at the point of discontinuity of 
the entire series. The convergence behavior of this latter expansion 
can thus be expected to be superior to the expansion in terms of ideal 
modes. 

The electric and magnetic fields of the imperfect asymmetric slab 
waveguide are expressed by the series expansions 

■ - £(**••,*•+ •"«,"> a) 

H = £ (c, (+ y (+) + **V*). (2) 
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The expansion coefficients c* +) and c F (-) are functions of the length 
coordinate z. The superscripts (+) and ( — ) indicate waves traveling 
in positive and negative z-direction. The sums in (1) and (2) are 
symbolic representations of a summation over guided modes plus an 
integration over the radiation modes of the continuum. 4 - 5 In order to 
simplify the notation, both sum and integral are indicated by the 
same symbol. In the integral, the summation index v is replaced by 
the continuous variable v, and the sum must be understood as the 
integral 

£-» f dv. (3) 

v J0 



The local normal modes E y and H y are solutions of the equations 

T#,(e, X 3fe, (±> ) + V, X 3C, (±) = za>e n 2 £, (±) (4) 

=Fi0,(e. X £» (±> ) + V, X S„ (±) = - iwnoKl^. (5) 

The upper and lower signs and superscripts belong together. The 
symbols appearing in these equations have the following meaning. 

j8, = propagation constant of mode v 
e« = unit vector in z-direction 
V< = transverse part of the operator V 
oj = radian frequency 
e = dielectric permittivity of the vacuum 
/i = magnetic susceptibility of the vacuum 
n = dielectric constant of the waveguide [n = n(x, y, z)~]. 

Substitution of the field expansions (1) and (2) into Maxwell's equa- 
tions and use of the orthogonality relations [see (9)] lead to the set 



x = f(z) 




CORE BOUNDARIES 
- x = - d+h(z) 

Fig. 1 — Sketch of the asymmetric slab waveguide with distorted core boundaries. 



(8) 
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of coupled wave equations 2 ,3,4 

j (+) 

dC " -a (+) , v- /^ (+ - +) (+) i v i+ '- ) (_ \ /a\ 
= — tfi^p,, + 2- (Km c " + &i» c " ) ( 6 ) 

dz , 

j (_) 

= i/3„c„ - L (#*.» c, + Kv, c, ). (7) 

dz , 

The coupling coefficients have the form 2,3 

/ (+ ). 3k, (+ \ I 

The asterisk indicates complex conjugation. The superscript p stands 
for (+) or ( — ) while the factor p assumes the values +1 and — 1. 
P is a normalization parameter which is related to the power carried 
by the modes via the relation 

- f (£„ (+) X H^Ve.dz = Pfi„,». (9) 

2 J—aa 

The symbol b vv > indicates the Dirac delta function if both v and v' 
represent continuous variables, it represents the Kronecker delta 
symbol if both v and v' are discrete labels, and it is zero if one subscript 
belongs to discrete modes while the other indicates a mode of the 
continuum. 

The coupling coefficients (8) are not very easy to evaluate since they 
are expressed in terms of derivatives of the mode functions. A. W. 
Snyder 3 has shown that the coupling coefficients can be transformed 
to the following more useful form. (/3„ (-) = — /3, (+) ) 

<±,, m ^ /■- a-^jB. cpi (10) 

4P(0<±) _ pM) ;_„ dz 

The coupled wave equations (6) and (7), with the coupling coefficients 
(10), provide an exact description of imperfect dielectric waveguides. 
The one-dimensional integral in (10) constitutes a specialization to a 
two-dimensional problem in view of our interest in the asymmetric 
slab waveguide. By extending the integration over the cross-sectional 
x, y plane, the general coupling coefficients are obtained. 
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For our purpose, it is advantageous to derive approximate coupling 
coefficients for asymmetric slab waveguides with discontinuous index 
distributions. We use the fact that the normal component E x , of the 
electric field obeys the relation 

2 2 

>iiE x 'i = n 2 E x > 2 . (11) 

It is shown in Fig. 1 that Wi and n 2 are the values of the refractive 
index at either side of the interface. In order to derive the desired 
expression for the coupling coefficient, we assume that the discon- 
tinuous index distribution is smoothed out (in an arbitrary way) into 
a continuous distribution. We assume that the wavelength of the 
radiation is very much larger than the region over which the refractive 
index varies continuously and write 

2 

E x . = -E«i. (12) 

n 2 

We show in the appendix how the integral in (10) can be evaluated and 
obtain the result 

2 

(±.p) we <» \, 2 idfVni (±) * ( P ) (±) * ( P ) 



K uv ' = { (wi — n 2 ) — — S.x 8 VX + &uv 8 



2 

2 



+ S uz 8 VZ — (»i — w 3)— —&ux 8 

Jx=/(0 dz \_n 2 



(p) 



(±)» (p) (±)» (?) 



lx d+h(z) 



■ (13) 



The index distribution can now again be considered as discontinuous. 
The functions /(z) and h (z) describe the deformation of the upper and 
lower side of the core boundary (see Fig. 1). The field components are 
taken inside the core region at the core boundary. The refractive index 
of the core is ni while n 2 and n 3 are the indices above and below the 
core region. The electric field components are related in the following 
way. 



si:' - 


Gyx 


8yy = 


8yy 


8y t = 


- £' 



(14) 
The approximation involved in the coupling coefficient (13) consists 
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in using the x component E, x and z component E vx of the local normal 
modes instead of their normal and tangential components with 
respect to the interface. The approximation is valid provided that 

df dh 

— «1 and — «1. (15) 

dz dz 

For many practical applications it is sufficient to use an approximate 
solution of the coupled wave eqs. (6) and (7). In particular, for the 
calculation of the radiation loss coefficient, we use the approximate 
solution of (6) (p = p) 

c p (z) = ± c v e-*"' * / K p , (u) 

•exp l-i f (pl'\v) - p ±y )dv\ du. (16) 

The coefficient ci p) is the amplitude of the guided mode, the losses of 
which we want to calculate, taken at z = 0. The propagation constant 
y is a function of z since it belongs to a guided mode in a non-uniform 
waveguide, /3 P is independent of z since it belongs to a radiation mode. 
The relative power loss AP V /P> that mode v suffers in traveling from 
z = to z = L is given in Refs. 5 and 6. 

AP r l r nik 



Py IC^I 2 



f m ' \c?\L)\Up. (17) 

JO 



The sum in front of the integral sign indicates that we must add up the 
contributions of forward and backward traveling modes as well as 
the contributions from the various kinds of radiation modes that will 
be discussed in the next section. The integral extends over the range 
of p values that belongs to propagating (non-evanescent) radiation 
modes. We show below that the functional form of the radiation 
modes is not the same over the entire integration range. 

III. MODES OF THE ASYMMETRIC SLAB WAVEGUIDE 

We consider only the special case in which there is no field variation 
and no waveguide distortion in y direction. This fact is symbolically 
expressed by the equation 

~ = 0. (18) 

by 



ASYMMETRIC SLAB WAVEGUIDES 69 

With the restriction (18), the fields of the slab waveguide can be 
classified as either TE or TM fields. 8 The TE fields have only the 
following non-vanishing field components 

E v , H x , H z . (19) 

The TM fields have the non-vanishing field components 

H v , E x , E,. (20) 

It is assumed that the refractive indices of the waveguide are ordered 
in the following way 

ni > n 2 ^ n 3 . (21) 

IV. GUIDED TE MODES 

The x and z components of the magnetic field follow from the E y 
component by differentiation 

i dEy 

H x = (22) 

uifio dz 

i dEy 

H z = -• (23) 

ojfio Bx 

The guided TE modes of the asymmetric slab waveguide are obtained 
as follows (the factor exp[i(o)t — /3z)] is always suppressed) : 

& v = A g e-* x for ^ j < co (24) 

Sy = Agl cos kx sin kx 1 for — d ^ x ^ (25) 

Sy = A ( cos Kd + - sin «d Je 9(x+d) for -*><x^-d. (26) 

The parameters appearing in these field expressions are defined by the 
equations : 

7 2 = 02 _ n \ k 2 (27) 

2 = /3 2 - n\k 2 (28) 

K 2 = n \w - /3 2 (29) 

k 2 = oj 2 6oMo. (30) 
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The propagation constant is determined from the eigenvalue equation 

tan /ca = (31) 

k 2 - yd 

The normalization of the mode is obtained by expressing the amplitude 
coefficient A g in terms of the power P carried by the mode. 

Al = ^f (32) 

|/3|(d + - + -V/c 2 + 7 2 ) 



The mode labels v, that were used in the coupled wave equations and 
the field expansions, are suppressed. The modes are labeled according 
to the solutions of the eigenvalue equation (31). 

V. TE RADIATION MODES 

The propagation constants of the radiation modes do not have a 
discrete set of values. However, the asymmetric slab waveguide has 
different types of radiation modes. In the range 

n 3 k g ^ n 2 k (33) 

we find only one type of radiation mode, the fields of which decay 
exponentially into the region three with refractive index n 3 , but are 
standing waves in the space above the waveguide with refractive 
index n 2 . We can visualize the physical mechanism for exciting these 
modes by assuming that a source at infinity sends a plane wave that 
impinges on the core of the slab waveguide under an angle that is 
given by 

cosa = — • (34) 

rtik 

The incident plane wave penetrates into the core region but is totally 
internally reflected if the angle a stays in the range given by (33). 
This total internal reflection occurs because we assumed that n 3 < n 2 . 
In the space above the core we find a reflected wave added to the 
incident wave supplied by the external source. This explains the occur- 
rence of standing waves in this region. It is not possible to find solutions 
of Maxwell's equations satisfying the boundary conditions which have 
only traveling waves outside of the core region. 

In the range of propagation constants given by (33) we find the 
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following expression for the field of the radiation modes 



a 



& y = A r cos px + -B r sin px for ^ x < « (35) 

P 

8 y = A r cos ax + i?r sin o-z for —d^x^O (36) 

<g„ = (A r cos ad - £ r sin ad)e-*' A(I+d) for - oo < x ^ -d. (37) 

H x and H, are obtained from E y with the help of (22) and (23). The 
constant B is related to the constant A by the equation 

A — ia tan ad 

Br = 77 7T- Ar (38) 

A tan era + ia 
and the parameters appearing in these equations are defined as follows 

ff 2 = nltf _ £2 (39) 



p 2 = n 2 k 2 - 2 (40) 

A 2 = n 3 k 2 - /3 2 . (41) 

Note that A is imaginary for ft values in the range given by (33). 

It is convenient to identify the parameter p with the mode label v 
to label the radiation modes. We thus use the identity 

*w = *(p ~ p) (42) 

in (9) and find for the amplitude coefficient A the relation 

/ A • V 

p 2 ( a cos ad H sin ad ] 

Al = ^ i * i (43) 

ir 1/8 1 / A V / , ^ A 2 

p 2 ( o- cos <rd H — - sin crd J -\- a 2 la sin ad 7 cos ad J 

With |8 in the range (33) we find that p is confined to the region 

^ p ^ (n\ - nl)*A. (44) 

Next we proceed to list the radiation modes that belong to propagation 
constants in the range 



0^/3^ n 3 k (3 real 
and 

0^ |/3 1 < °° /3 imaginary 



(45) 
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The corresponding range for p is given by 

(nl - n\)*k < p <°o. (46) 

For real values of /?, these modes propagate along the z axis while they 
are evanescent waves in z direction for imaginary values of /3. It is 
again possible to visualize the modes in the range (4G) as being excited 
by a source outside of the waveguide core located at infinity. This 
source sends a plane wave toward the slab whose angle of propagation 
with respect to the z axis is given by (34). However, there is now no 
longer total internal reflection at the lower boundary of the core so 
that we obtain an incident and reflected (in x-direction) wave in the 
space above as well as inside the core. Below the core there is a trans- 
mitted propagating wave. However, we may now assume with equal 
justification that a second source sends a plane wave in the direction 
of the core from below. If both sources are turned on simultaneously, 
we obtain standing waves (in x-direction) below as well as above the 
waveguide core. The exact form of the radiation field depends on the 
relative phases between the two sources. It is thus not surprising that 
there should be an infinite number of ways in which orthogonal sets 
of radiation modes can be constructed. 

We list only the E y components of the modes and obtain the H x and 
H, components by differentiation from (22) and (23). (i = 1, 2) 

S„ = C r ( cos P x + -Fi sinpxj ^ x < <*> (47) 

S v = C r (cos ax + Fi sin ax) -d ^ x g (48) 

&v = Cr\ (cos ad — Fi sin ad) cos A (x + d) 



-\ (sin ad + Fi cos ad) sin A (x + d) \ — °o < x ^ — d. (49) 

A J 

The parameters a, p and A are given by (39), (40) and (41), A is now 
a real constant. Whereas the amplitude coefficient B T in (35) through 
(37) was related to A T by the boundary conditions, we now face the 
situation where the amplitude coefficient Fi remains arbitrary. Equa- 
tions (47) through (49) satisfy Maxwells' equations as well as the 
boundary conditions without any further restriction having to be 
imposed on the coefficient Fi. This freedom of choice is related to the 
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arbitrary amplitude and phase of the two plane wave sources of the 
radiation mode mentioned above. We must choose Fi in such a way 
that two radiation modes with the same value of the propagation 
constant, but different values of Fi become mutually orthogonal. But 
even this requirement does not specify the possible values of Fi 
uniquely. We are thus free to choose F values according to our own 
convenience. Of the infinitely many possibilities, we choose the F x 
coefficients such that in the limit n 2 = n 3 , even and odd radiation modes 
result. In the asymmetric slab waveguide no even or odd modes exist. 
But the guided modes become either even or odd as n 2 approaches 
n 3- We obtain the same symmetries for the radiation modes by a 
suitable choice of the Fi coefficients. 

Fi ' 2 = T* Att • 9 J (ff2 ~ ^ C ° S ^ + ~ (ff2 " ^ 

(a 2 — A 2 ) sin 2ad I p 

± (a 2 - A 2 ) 2 4- 2— {a 2 - A 2 )(<r 2 - p 2 ) COS 2ad 
P 

A 2 1*1 

+ -(* 2 - P 2 ) 2 • (50) 
P 2 J I 

The + sign (— sign) belongs to the odd (even) mode in the limit 
n% = ft3, A = p. We identify the mode label v again with the transverse 
propagation constant p, and obtain from the normalization condition 
(9) and (42) the relation between the amplitude coefficient C T and the 
power P 

2 4a)/xoPf A 

C T = — ■ —(cos ad — Fi sin ad) 2 

r\p\ [_p 



4- —(sin ad 4- Fi cos ad) 2 + 1 + — F, 2 
Ap p 2 . 



(51) 



We have thus obtained two distinct sets of radiation modes, the 
propagation constants of which lie in the range (45) . The two sets are 
distinguished by the labels i — 1 and i = 2. F\ and F z follow from (50) 
if we use both signs of the square root expression. It is noteworthy that 
the following relation holds. 

F t F 2 = -1. (52) 

This listing contains the complete set of TE guided and radiation 
modes consistent with the restriction (18). All modes are mutually 
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orthogonal to each other. We have concentrated on the forward 
traveling modes. The backward traveling modes are obtained simply 
by changing the sign of 0, (/3 (_) = — /3 (+) ). 

VI. GUIDED TM MODES 

The TM modes are very similar to the TE modes except that the 
roles of the field components are interchanged. We now list only the 
H v component and obtain the E x and E, components of the field by 
differentiation 



(53) 

(54) 

(55) 



i BH V 

F 




nrwea dz 




-i dH v 

F 




rruea dx 




3C„ = DgBTT 


^ X < co 


3C„ = Dgl cos kx 


2 

sin kx 


\ 


n 2 k 

2 



) 



-d ^ x ^ (56) 



2 

(ni t \ 

cos Kd H sin nd \e^ x+d) 
n 2 k / 



- oo < x ^ -d. (57) 



The parameters k, y and are defined by (27) through (29). The 
eigenvalue equation is 

tzn K d = nUnly + nl6) (58) 

n*nV - n 4 y6 

2 3 1 

and the amplitude coefficient is given by 
4coe P 



D„ = 



1/31 

i 4 bv r njnz k 2 + y 2 , wina k 2 + 2 1 
(n 2 /c 2 + tixy 2 ) d + + — 



Win2»c 4 , % 

X- — — (59) 

2 2 ... 22 



1 3 1 
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VII. TM RADIATION MODES 

The TM radiation modes must again be split up into two ranges. 
In the range ^ p ^ (nj — n|)*fc, the fields have the form 

2 

ni a 

3C y = D r cos px H G> sin pz ^ a; < « (60) 

n 2 p 
i 

3C„ = D, cos ax + G T sin arc —d^x^O (61) 

3C„ = (D r cos ad - G T sin <rd)e-<* <*+<*> - oo g x ^ -d. (62) 

The boundary conditions require that G> is related to D T in the following 
way: 

2 2 

7iiA cos ad — inz<r sin ad 

G r : : D r . (63) 

w 2 A sin ad + in 2 a cos ad 

1 3 

The parameters a, p and A are defined by (39) through (41), — i'A is 
a real positive quantity. The amplitude coefficient D T is related to 
the power coefficient P 

2 

2 WlzUtoP 1 

D r = — (64) 

r\S\ 4 2 x ' 

n 4 p 2 D 2 

In the range (n| — n|)*Ai ^ p ^ » the radiation modes have the form 
(* = 1, 2) 

2 

3C„ = £ r ( cos px -\ Ri sin px ) ^ x < w (65) 



cos px H ft, sin pa; J 

njp / 



0C„ = *S r (cos ax + fti sin ax) -d ^ x ^ (66) 



3Ci/ — Or 



(cos ad — Ri sin ad) cos A(a; + d) 



2 

Wi * . 

H (sin o-d + fti cos trd) sin A (a; + d) 

n 2 A 
i 



- oo < x g -d. (67) 
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The coefficients Ri are again arbitrary. Our choice 

1 [44 

Ri o = { (n 3 a 2 — ??iA 2 ) cos 2ad 

(nV - n 4 A 2 ) sin 2ad 

S 3 1 

n-z A « 4 f 4 4 n 3 A 2 4 4 

+ (n 2 a 2 - nip 2 ) ± (rc 3 <x 2 - niA 2 ) 2 + (n 2 o 2 - n lP 2 ) 2 

n 2 p L n l p 2 

2 
7*3 A 4 4 4 4 ~|'] 

+ 2 (?i 3 a 2 - «! A 2 ) (n a «r 2 - nip 2 ) cos 2ad (68) 

nj P J J 

causes the modes with index i = 1 to be orthogonal to the modes with 
index i = 2 and, in addition, assures that these modes become even 
and odd in the limit n 2 = n 3 . [The + sign (— sign) belongs to the 
odd (even) mode.] The normalization is given by 



2 4we P [ 1 ii 2 ff 2 2 n 3 a 2 

S r = \ 1 R< -\ (sin ad + Ri cos ad) 2 

t \p\ In 2 n 4 p 2 n 4 P A 



2 1 



1 A l" 1 

-\ (cos ad - Ri sin ad) 2 \ . (69) 

n l p j 

All amplitude coefficients for the TE and TM modes were taken to 
be real quantities. This assumption does not lead to a loss of generality 
since the necessary phase factors are incorporated in the expansion 
coefficients c,. 

VIII. COUPLING COEFFICIENTS 

With the help of the expressions for the normal modes and the 
coupling coefficients (10), any problem of the asymmetric slab wave- 
guide with arbitrary irregularities of its refractive index distribution 
can be solved, provided that the restriction (18) is imposed. Problems 
caused by core boundary irregularities or by gentle tapers can be solved 
with the help of the coupling coefficients (13). For convenience, a few 
coupling coefficients will be worked out explicitly. 

As long as the restriction (18) applies, TE modes do not couple to 
TM modes. All coupling coefficients between TE and TM modes 
vanish. We restrict the discussion to listing the coupling coefficients 
between guided TE modes, guided TM modes, and to coupling from a 
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guided TE or TM mode to its respective radiation modes for the case 
of core boundary irregularities. 

The coupling coefficients between two guided TE modes can be 
obtained from (13) and (25). (p = + or — ) 



K 



(±.p) 



[df sin k 
t ~ ~r- — 
dz sin k 



sin K u d sin K„d dh" 



4 sin K v d I dz_ 



(0, - 



(p) 



)[ i w, i (a + - + -Yd + - + ~ yi 



(70) 



The eigenvalue equation (31) was used to express (70) in this simple 
form. This coupling coefficient (and all others to be listed below) holds 
for the special case /(z) <3C t/k, h(z)<£ tt/k with k of (29). Instead of 
using the values of the field at x = f(z) and x = —d-\-h(z), the field 
values at :r = and a; = —d were used. In order to see the agreement of 
this coupling coefficient with the coupling coefficient for the symmetric 
case [eq. (7) of Ref. 7], it is necessary to note that the core thickness 
d of this paper corresponds to 2d of Ref. 7. In addition, we need to 
keep in mind that only the Fourier components of f(z) and h(z) 
with spatial frequency /3^ ±) — & p) contribute to coupling between 
modes /i and v. The derivatives appearing in (70) are thus equivalent 
to the products i(fi^ - /3< p) ) /(z) and itf^ - pi p) )h(z). Keeping 
these remarks in mind, complete agreement of (70) with (7) of Ref. 7 
is obtained for the special case n-i = n 3 , y = 6. 

The coupling coefficient for coupling between a guided TE mode v 
and a TE radiation mode p follows similarly from (13), (25), and (36) 
for radiation in the range ^ p S {n\ — n|)*A; 



K 



(±,p) 



•iii A 

(wi — n 2 ) } fc«cp( o- cos ad + -7 sin ad 



(0(±) _ 0<p)) 



2 2 
r df sin K,d (n\ — n?\ * 



dh 



Idz I sin n,d \ \n 2 — n 2 / A , dz 

1 2 a cos ad H sin ad 

i 

• I ir I /3 p /3„ 1 1 d H 1 }| p 2 ( o- cos ad -\ sin ad 



7, 0, 



(• 

( <r sii 



+ o- 2 ( o- sin ad -cos ad) \\ . (71) 



78 THE BELL SYSTEM TECHNICAL JOURNAL, JANUARY 1973 

The coupling coefficient of the guided TE mode v to the TE radiation 
mode p in the range {n\ — n\)*k ^ p < <» is given by 

2 2 

(±.p> (»i— «»)***» 

Kj pv = 



(0(±> _ 0<p)) 



7r|(8 p /3,|(d+ — + -) 
\ 7* ft/ 



2 2 

d/ sin K v d (n\—n%\l . d/i 

(cos ad — F j sin ad)- 



dz I sin <c„d I \n 2 — n 2 / dz 



TA <r 2 er 2 2 -|» 

—(cos ffd-Fj sin «2)H (sin <rd+Fj cos crd) 2 +H — F> 

Lp pA P 2 J 



(72) 



The factors Fi and /^ are obtained from (50). The radiation modes do 
not all propagate along the z-axis. Propagating modes are confined to 
the range ^ p ^ n 2 /c. 

The reader should not be startled by the fact that the coupling 
coefficients (70) have the dimension m~ l while the coupling coefficients 
(71) and (72) have the dimension mr*. The different dimensions are 
attributable to the fact that the coupling coefficients between guided 
modes occur under a summation sign while the coupling coefficients 
that describe coupling to radiation modes occur under an integral sign. 
The integration is performed with respect to p, the dimension of which 
is mr 1 . The product c„K pv dp has the dimension m~ l in agreement with 
the dimension of the coupling coefficients for guided modes. 1 " 

Finally, we list the coupling coefficients for the TM modes. Coupling 
between guided TM modes is described by the coupling coefficient 
(p = + or -) 

2 2 

_^(±.p) (»i — n2)D gv D 0ll ( 2 (±) (p) 2 df 

4P(/3<±> - /3<»>)co6 7i 2 n 4 l v dz 



1 2 



J • J 2 2 r-( 4 2 I 4 \( 4 2 I 4 2 \-l* 

sin K v d sin K M a n\ — n 3 \{n^c p ■+■ Mi7 ,..Hn 2 K M -+- ?k7*J r 



[~(n 2 K, + ni7„) (n 2 K M + Wi7m)"| ! 
L (n 4 « 2 + n 4 2 )(n 4 K 2 + n 4 2 ) J 



I sin K„d sin K„d| n 2 - ra 2 L (n 4 « 2 + n*fl")(nV + n 4 . 

1 M ' 1 ! J» 1 / 3 ft 1 M 

5 (4-1 («1 2 CS/i. I 

•(«rfj ft +n,W^T • (73) 
azj 

The coupling coefficients for the TM modes are far more complicated 

t Note that /* |c p l 2 d p is dimensionless so that the dimension of c„ is w*. 
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than the corresponding coefficients for the TE modes. To simplify the 
notation, we did not substitute expression (59) for the mode amplitude 
into (73). 

The coefficient for coupling from a guided TAX mode to a TM 
radiation mode in the range fS p ^ {n\ — n\) h k is 

2 2 

<±,p) (ni - n 2 )Dg t D ri , |/_<±>_ci.) G rp \df 

K DV = — — — - UP, P, - °y» ^r~ 



4P(/3 (±) - ^ (p) )«eo« 8 n 2 l\ ' D Tp J dz 

p v 12 

2 2 4 2 4 2 

sin K v d tti — n 3 fn 2 Ky + ni7„~|* I" ( ± ) < p > / G> 



— nlrnlicl + ni7,"| ! [~ <±> < P ) / , G rp . \ 

\ \ B„ j8„ cos <rd sin crd 

|sinjc,d| n 2 - n 2 LnV + n 4 2 J L \ D rp J 

1 ' 1 2 3 v If 

/ G rp \~\dh] , x 

4- <rfl, ( sin ad H cos <rd J — ■ (74) 

The amplitude coefficients D„„ D Tp and G rp are obtained from (59), 
(63) and (64). 

The coefficient for coupling from a guided TM mode to a TM 
radiation mode in the range (n\ — ni) ifc ^ P ^ °° is given by 



(± lP ) (ni — n 2 )D gv Sr 

Ki pv = 



dz 



2 2 4 2 4 2 

sin K„d ni — nzVntft, ■+■ ni7»T 



I sin K„d I n 2 — 

1 ' 1 2 



4P( j g(±) _ ^p))we nV 

2 4 2 4 2 

W 3 rr?2K, + ?U7vl J I" (±) (p) . 

— /3 P j8 r (cos ad — Ri sin ad) 

n 2 U 4 k 2 + n 4 2 J L 

2 3 v I v 

Idh) 
+ a0„(sin <rd + R, cos ad) — ■ (75) 

The amplitude coefficients D gv , S rp , and Ri are obtained from (59), 
(68), and (69). The index i assumes the values 1 and 2, and corresponds 
to the two types of radiation modes that are distinguished by the + 
and — signs in (68). The superscripts, + and — , attached to the 
coupling coefficients are supposed to indicate whether the modes 
travel in + or — z-direction. 

It can be shown that the coupling coefficients derived in this paper 
specialize to the correct expressions 4 - 5 of the symmetric slab waveguide 
in the limit n 2 = n 3 , y = A, p = 0. f 

t Equations (9.5-26) and (9.5-27) of Ref. 5 must be divided by n|, eq. (9.5-31) 
must, correspondingly, be divided by nj. 
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IX. CONCLUSIONS 

We have collected the formulas for the modes of the asymmetric 
slab waveguide and have used this information to derive the coupling 
coefficients between guided modes as well as between guided and 
radiation modes for the case of very slight core boundary imperfections. 
Also presented is the general theory of coupled modes of dielectric 
waveguides and the general formulas for the coupling coefficients. 
The theory collected in this paper is useful for the description of mode 
conversion and radiation phenomena. Phenomena such as the grating 
coupler and the interaction of acoustic waves and guided light waves 
can readily be treated with the theory presented here. For an applica- 
tion to statistical irregularities of the core boundary, the reader is 
advised to consult Refs. 5, 6, and 7. 

APPENDIX 

Evaluation of the Integral {12) 

We consider the index distribution of the slab waveguide as being 
smoothed out to avoid the discontinuity at the core boundary. It is 
assumed that the index varies only in a direction perpendicular to the 
core boundary. We use a coordinate system x', z', the axes of which are 
perpendicular, and parallel to the tangent at a particular point of the 
core boundary as shown in Fig. 2. In this coordinate system, we assume 
that the refractive index is of the form 



n 2 = F(x'). 



(76) 



For values of x' ^ we have F = n\\ for values x' > £ we have 
F = n\. At the end of our discussion, we let £ — > 0. The scalar product 




TANGENT OF 

CORE BOUNDARY 



Fig. 2 — Sketch of the coordinate systems used for the evaluation of integral (10). 
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of the electric field vectors can be expressed as 

Sj-fi, = C"^*< +<**»«■ (77) 

The coordinate £ indicates a direction parallel to the core boundary. 
S, t is continuous at the core boundary and can be taken out of the 
integral. S» X ', on the other hand, is discontinuous. We express it in 
terms of the field just inside of the core and, using (12), write 

2 

&„• - ^r(«"')«-o. (78) 

F(x') 
The scalar product can thus be written in the form 

4 
£ M - £„ = — (£ M i'<§,x')x'=0 + SptSyt. (79) 

F 2 (x) 

We thus have to deal with two different integrals. We first consider 
the integral 

B ^,^r^ d , (so 



(81) 



dz dz J o dx' 

We obtain from Fig. 2 the relations 

dx' 

— = sin a 
dz 

and 

dx = —. (82) 

cos a 

The integral can thus be evaluated 

r°° dF(x') 

/ — — dx = [F({) - F(0)] tan a . (83) 

Jo dz 

At the upper core coundary, we have tana = df(z)/dz, and at the 
lower core boundary we have tana = dh(z)/dz. Taking both core 
boundaries into account, we find with the help of (76), 

/" dn 2 . 2 2 df * 

— SatSvtdx = — (ni — n 2 ) — (5 Ml 5 w )*=/ 
.«, dz dz 

2 2 rfft • 
+ (ni — n 3 ) — (8m8 v t)x — d+h- (84) 
dz 
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The integral associated with the normal field components is essentially 
of the form 

1 W J u J" 1 d *W J / 
ax = (tan a) I dx 



F 2 (x') dz 



F(£) - F(0) 

= (tan a) — — • (85) 

F(0)F(S) 

The integral containing the normal field components results in 



/. 



00 dn 2 . 22 Mi . df 

Sr X >8 rX 'dx = — (ni — n>) — (©„«'&»')»-/ — 

_„ dz n* dz 

2 



2 2 «i * dh 

+ (wi — n 3 ) — (&?*'$,*')* — d+h —■ (86) 
n 2 dz 

3 



In (13) we replaced x' with x and £ with z. These approximations are 
valid provided that the inequalities (15) apply. The error is of second 
order in the derivatives of /(z) or h(z). 
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